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1. Introduction 



S 
> 

K> \ We consider estimation of the p-dimensional location parameter of a spherically symmetric 

C^ ■ distribution. Specifically, let X = (Xi, . . . ,Xp)' have a density function f{\\x — 9\\) and 

consider estimation of 6 with a general quadratic loss function Lq{9, 5) = {5—6yQ{5—0) = 
\\5 — 0||q for a positive definite matrix Q. The usual miniinax est imator X, which is 
genera lized Bayes , is inadmissible for p > 3 as shown in ISteinI ( 19561 ) in the normal case, 



and in iBrownl ((l966) under more general situation, respectively. In the decision-theoretic 



point of view, we are interested in proposing admissible estimators dominating X, that 
is, minimax admissible estimators. Note that the dominance over X means minimaxity in 
our setting because X is minimax with a constant risk. Note also that our results hold 
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for the location vector of elliptically-contoured distributions, because we are considering 
a general quadratic loss function with arbitrary positive definite Q. 

In the n o rmal case, there already exists a broad class of admissible minimax estimators. 
BaranchikI ( 19701 ) gave a sufficient condition for minimaxity of a shrinkage estimator of 
the form 

M^) = (1 - mx\\)\\x\\-')x. (1.1) 

StrawdermanI ( 197lh found a subclass of pro per Bayes esti mators of the form (|l.lj) satisfy- 
ing the sufficient conditions for minimaxity. iBrownl ( 197l|) gave a very powerful sufficient 
conditi on for a dmiss ib ility of generalize d Baye s est imators. Us i ng Brown's (1971) con- 
dition, iBergen (|l97a), JFourdrinier et al.l (|l99^) and iMaruvamal (|l998l . |200J) enlarged a 
class of admissible minimax estimators which are of the Strawderman type and gener- 
alized Bayes. For a subclass of scale mixtures of multivariate normal distributions which 
inc ludes multivariate -t distribution, some proper Bayes minimax estimators were proposed 
by IMaruvamal (2003) by using Strawderman's (1971) techniques. 

However, for general spherically symmetric distributions, no minimax admissible esti- 
mators of the location vectors have been derived, although for the mini maxity. variou s 
Baranchik-type sufficient conditions of the estimator (11.11) were given by Bergei (|l975|), 
Brandwein and StrawdermanI ( 19781 . Il99lj ) and iBockl f|198fii) . The main reason is the lack 
of a standard class of generalized or proper Bayes estimators of the form (jl.lj) like the 
Strawderman type in the normal case, which allows an easy check of the minimaxity condi- 
tion. Furthermore n o sufficient condition for admissib ility of generalized Bayes estimators 
has been derived. In lMaruvama and Takemural (|200fJ ). we provided satisfactory solutions 
to these problems. In this pape r we will weaken some regularity conditions assumed in 
Maruvama and Takemural (|200.^ . 

In Section 2, we give preliminary results including the properties of regularly varying 
functions and asymptotic behaviors of expected values when \\9\\ is sufficiently large. The 
former is useful for constructing a very convenient sequence of proper densities approaching 
an improper density g{9), which is required in applying the method of Blyth (1951). 

In Section 3, we will present a powerful sufficient condition for admissibility of gener- 
alized Bayes estimator and in particular show that the generalized Bayes estimators with 
respect to the harmonic prior g{9) = \\0\\'^^'^ and with respect to a prior with a slightly 
heavier tail 



9iO) 



'log( 



+ c 



c> 1, 



(1.2) 



are admissible under mild regularity conditions on /. 

In Section 4, we show that the generalized Bayes estimator with respect to the harmonic 
prior is written as 

j,Hp-'F{t\\x\\)dt' 



1 



X, 



1 I-, (1-3) 

V J,'tP-^F{t\\X\\)dtJ 

where F(u) = J^ sf{s)ds. This form is simple enough to check various sufficient conditions 
for minimaxity and we demonstrate that ()1.3() is minimax for some /. We believe that ()1.3() 
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is niinimax for a broad subclass of spherically symmetric distributions. Notice that the 
generalized Bayes estimators with respect to priors except ||0|P~^ do not have such simple 
forms as far as we know. 

Ou r proof of admissibil i ty is mainly based on the techniques of iBrown and Hwand 



( 19821 ) . iBrown and Hwand ( 19821 ) considered the problem of estimating the natural mean 



vector of an exponential family under a quadratic loss function. Note that the intersection 
of their setting and our setting is the multivariate normal case. Their sufficient condition 
for admissibility in the normal case does not however permit g(9) to diverge to infinity 
around the origin l ike ||0|P~^, while it permit s q(9) < ||0||" with a < 2 — p for sufficiently 
large \\9\\. Prior to lBrown and Hwangt (1982 ). lBrownl (|l97ll ) considered the estimation in 



the multivariate normal case and gave a powerful sufficient condition for minimaxity which 
are satisfied by the harmonic prior and l|1.2|) , but his proof was based on r nanv advanced 
mathematics. Our mathematical tool is much more familiar to the readers. iBrownl ( 197lh 



also gave sufficient condition for inadmissibility. By using it we see that the generalized 
Bayes estimator with respect to ||^|p~^log (||^|| + 2) is inadmissible. Hence our sufficient 
condition for ad missibility should be very tight and close to being a necessary condition. 
Brownl (|l979l ) considered a more general problem than ours: estimation of 6 for a general 



density p{x — 9) and a general loss function W{6 — 9). He conjectured that the prior 
g{9) ~ ll^ll" with a <2 — p for sufficiently large \\0\\ le ads to admiss ibilitVi regardless of 
the density p and the loss W . Hence our results support JBrownl ( 1979|)'s conjecture for the 



case of elliptically-contoured family and a general quadratic loss function. 

Finally we notice that the most important key for our proof for admissibility is the 
constructi on of a very convenient se quence hiiO) for approximating g{9) by g{9)h]{9) for 
< 7 < 2. JBrown and Hwang! (|l982r used 7 = 2 and 



hi{9) 



1 \\e\\ < 1 

l-log||6i||/logi l<||6l||<i 
ll^ll>i. 



This hi{9) is not differentiable at \\9\\ = 1 and truncated at ||^|| = i, which makes handling 
and extension difficult for our purposes. Our hi{9) given in Section 2 is smoother and not 
truncated. Furthermore our 7 is flexible whereas 7 = 2 in lBrown and Hwand ( 19821 ). By 



such a flexible 7, we can adjust the rate of convergence of g{9)h]{9) so that it is just 
enough to be proper. We will see that choosing as small 7 as possible is important in 
the main theorem, Theorem 3.1. We believe that our smooth function hi{9) and an idea 
of flexible 7 are very useful for showing admissibility of generalized Bayes estimators in 
various problems. 

2. Preliminaries 

In this section we prepare a sequence of proper densities using the theory of regularly 
varying functions and give some results on asymptotic behaviors of expected values when 
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the location parameter diverges to infinity. For the theory of regularly varyin g and slowly 



varyin g functions the readers are referred to lGeluk and de Haan (1987 ) and Bingham et al 
(1983). 



2. 1 . Regularly varying functions 

A Lebesgue measurable function / : i?+ —i- R which is eventually positive is called regularly 
varying if for some a (z R 

lim ^^ = r, Vt > 0. (2.1) 

We sometimes use the notation / € RVq. The number a in the above is called the index 
of regular variation. A function satisfying (|2.1|) with q = is called slowly varying. 

Let P : -R"*" -^ R^ be positive, continuously differentiable, monotone decreasing, inte- 
grable (i.e. J^ f3{r)dr < oo) and regular varying with index —1. A typical (3{7]) is 



1 1 "-^ 



rj + c {Log„(7? + c)}2 f-J-^ Logi(r/ + c) ' 
where n is a positive integer, Logo(?7 + c) = 1, 

Logj(7? + c) = log log • • • log(r/ + c), i > 1, 
^ ^ ^ 

i 

and c is chosen such that Log„(c) > 0. Note that for ()2.2|) 

1 



n T^^Y.^.v (2-2) 



ri Log„(?7 + c) 

The following results for (3{rj) satisfying the above assumptions are known from the 
theory of regularly varying functions. 

Lemma 2.1. 1. J^ (3{r)dr G RVq and I3'{r) € RV_2. 

2. limr,^^ r]f3{7]) / J^ p{r)dr = 0, lim^_.oo r?/?'(r?)//3(r/) = -1. 

We now define functions Hi{r]), i = 1,2,..., based on /3(r) by 

H^in) = '" .oo ., ' • (2.3) 

These functions are very useful for constructing a sequence of proper prior densities ap- 
proaching the target improper density in the next section. The properties of Hi are given 
in the following theorem. 

Theorem 2.1. i. < Hi{7]) < H2{r]) < • • • < 1. For any fixed t], limj_^oo Hi{'q) = 1. 
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2. For any fixed i, lim^^oo /„°° I3{r)dr f3{rj)~^ Hi{rj) = i and hence Hi(j]) G RV_i. 

3. For any fixed rj, limj_>oo H[{ri) = 0. 

I \H[{ri)\ < 2p{ri)/J^P{r)dr for all r? > 0. 

5. For any e > 0, there exists tjq such that —1 — e < r]H'-{r]) / Hi{r]) < for all r] > tjq 
and for all i. 

Proof. It is obvious that < Hi{rf) < 1 and Hi{rj) is increasing in i. For fixed r/, Hi{r]) ] 1 
by the monotone convergence theorem. 

By integration by parts, the numerator of Hi[rf) is written as 



oi'n-r)/i 



P{r)dr = i(3{rO + i j e 



{ri-r)/i(Dl 



f3'{r)dr. 



Therefore 



H^{ll) 



Piv) 



'I.^mdr 



J^°°e(^-'^)/^/3'(r)dr 



(2.4) 



(2.5) 



(^31) divided by (ESJ is 



H„ 



J^e-'/'p'{r)dr 
iv) I,^ I3{r)dr ' ' f^^ e-/^P{r)dr ' 



Piv) 



+ i 



For fixed i, the second term of the above equation converges to as ry ^ oo by the 
L'Hospital theorem. Hi{r]) G RV_i because J°° (3{r)dr £ RVq and /3(r/) € RV_i. 
Using (|2.4|) again, differentiation of the numerator of Hi[rj) gives 



,{r]-r)/i 



(3{r)dr 



,(v-r)/i 



I3{r)dr — (5{rj) = j e 

V 



{r)~r)/inl 



I3'{r)dr. 



Therefore 



HiiTj) 



I3{v) XT e^''~"'^/'(^{r)dr J^ e^''-''^/'{-P'{r)}dr 



(/r /3(Orfr)2 



I^ Pir)dr 



(2.6) 



Note that —(3'{r) > by our assumption. Each term of the right hand side of 1)2. 6|) is 
nondecreasing in i and hence by the monotone convergence theorem 



hm H[{rj) 



P{v)!:^mdr !:^{-P'{r)}dr 
U^I5{r)drY f^P{r)dr 



0. 



Furthermore we have 
\H'M\ < 



P{r])J^e^''-'y'P{r)dr 



U;°P{r)dr) 



+ 



J°°e^^-'-y'{-p'{r)}dr 



rPir)dr 
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POO 

< 2(]{r])/ / I3{r)dr. 



Dividing ((T^ by (|T^ . we have 

vPiv) C e--/'{-r/3^(r)//3(r)}/?(r)dr 

^/;°/3(r)dr /^°°e-A/3(r)dr 

r](3{r]) —rj3'{r) 

By 2 of Lemma l2. II the right hand side converges to —1. This impUes that for any e > 
there exists r/o such that r]H-{ri)/Hi{ri) > —1 — e for ah r] > r]Q and for ah i. Finally we 
will prove that H[[r]) < for sufficiently large r] independent of i. By 2 of Lemma l2.H 



vI^P{r)dr f f3{^) y_ /?'(r,) f3{i^) 



and hence (3{rj)/ /°° [3{r)dr is eventually nonincreasing. Hence by redefining r/o if necessary, 
we can assume that f3{r])/ J°° (3{r)dr is monotone nonincreasing for r] > rjo- By integration 
by parts on the numerators of each term in (|2.6() , we have 

e-''/'l3{r)dr = e"''/* / f3{r)dr - r^ / e^"^/* / /3{s)ds \ dr 

OO /-OO 






and hence 






which is nonpositive for rj > ijq. D 

2.2. Asymptotic behavior of expectations 

In the next section, we need evaluation of an asymptotic behavior of expectation 

EApm= I p{e)f{\\e-x\\)de 

JRP 
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for sufficiently large ||x||, where a random vector has the density function f{\\9 — x\\). 
This is the expected value with respect to the posterior distribution. Interchanging the 
roles of X and 9, in this subsection, we consider the asymptotic behavior of expectation 

Ee[p{X)]= I p{x)f{\\x-e\\)dx 

for sufficiently large \\9\\, where a random vector X has the density function f{\\x — 9\\). 
We believe that this does not confuse the readers. 

We discuss some notations used in the following. In addition to the Euclidean norm 



\x 



P = xf + • • • + Xp, we consider the norm ||x||^ = df xf + • • • + dpXp. For convenience we 
assume, without loss of generality, that di > ■ ■ ■ > dp > 1. Under this assumption 

\\x\\ < \\x\\d < di\\x\\. (2.8) 

By introducing this norm our results hold for elliptically-contoured distributions. The 
gradient of p(x) is denoted by 

Vp(x) = (Ap(,),...,_^^(,)y. 

We also write Vjp{x) = {d/dxj)p{x). Finally we write Cp = 2TrP'^/T{p/2). 

Now we make the following regularity conditions on the density / and the function p. 

Fl There exist tq > 0, L > 0, and s > 1, such that rP'^^f{r) < L for all r > tq. 
Bl p{x) is written as p{x) = ^(||x||d), where g^r) is continuously differentiable in r > 0. 
B2 There exists ri > 1 and ti < t2 such that Q(r) > and ti < rQ'{r)/ Q{r) < t2 for all 
r >ri. 

Assumption B2 is, for instance, satisfied by 

f f^ ~P + ct cos t \ 

Q[r) = exp / dt , for a > 

\Jo t + 1 J 

where we easily see rQ'{r)/ Q{r) = {r/{r + 1)}(— p + acosr) and hence ti = —p — a and 
^2 = —p + a in B2. Since lirur-^oo rQ'{r)/ Q{r) exists for regularly varying g, we deal with 
a broader class of functions than the class of regularly varying functions. We will discuss 
more in Section 3. Note that 

vT (- \ 6 \\\x\\d) ( -a ,2 \i 

Vp(x) = —r—r (d^Xl , . . . , dpXp) 

\m\d 

and 

||Vp(x)|l = k^^MM(4x? + . . . + 4^2)1/2 < dl|^'(||x|I,)|. 

W^ II" 
The following lemma is useful. The proof based on the integration of (log^(r))' = 
Q'{r)/ Q{r) is easy and omitted. 
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Lemma 2.2. Under the assumption B2 

{z/yY' < Q{z)/Q{y) < {z/yY^ 
for any z > y > ri. Moreover 

limsup sup ^(z)/^(y) < max(a*\/3*^) 

j/^oo ay<z<f^y 

for any < a < 1 < /3. 

We now state the following theorem concerning the asymptotic behavior of Eg[p{X)] 
for large ||0||rf. 

Theorem 2.2. Assume Fl, Bl and B2. For a = or 1, and j = l,...,p, if s > 
max(l, — ti — a — p,t2 + a) and /q r'^^^^^\g{r)\dr < cx), then there exists e > (say 
e = min(l, s + ti + a + p)/^) such that 



\e—a 
\d 



Ee[X<^p{X)]-e'^p{e) <Cp{9) 



(2.9) 



for \\9\\d > 2di max(ro,ri). Moreover C depends on p (or g) only through ri, ti, ^2 o.nd 
{eiri)}-^IS'rP+--^\e{r)\dr. 

For simplicity, in the rest of the paper, we will write E^lX^ p(X)] ~ 9'jp{6) if, as in 
Theorem 12.21 there exists e > such that ()2.9|) is satisfied for sufficiently large \\0\\fi- 

Proof. Note that 

Ee[X'^p{X)] - e^p{9) = Ee[X^{p{X) - p{9))]. 

Fix < z^ < 1 (set 1/ = 1/2 finally). Define 

K = {x : \\x - 9\\ < u\\9\\} 

VI = {x:(l-z.)||0||<||x||<(l + z.)||0||} 

K" = {x:d^\l-u)\\9\\d<\\x\\d<di{l + u)\\9\\d}. 



By ((TBI) K C VI C V//. Then 

E[Xf{piX)-pi9))] 



\e—a 
Id 



< 



\e—a 
\d 



+ I ]\\xr\p{x)-pi9)\f{\\x-9\\)dx 



< Md sup 



X d 



\pix)-pi9)\fi\\x-9\\)dx 



\€—a 
Id 



p{9) \\xrf{\\x-9\\)dx + 
v,9 



\e—a 
\d 



y,9 



\xr\p{x)\f{\x-9\)dx 



h + h + h. (say) 



(2.10) 
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Consider the first integral /i. sup^^y,, {\\x\\d/\\d\\dT ^ "^1(1 + ^)"- If s > 1, then mi 
Irp in ~ ^11/(11^ ~ 6\\)dx is finite. Therefore for \\9\\d > di(l — i')^^ri we have 

prJ \p{x)-pmf{\\x-e\\)dx 

= m\d I \{x - eyvp{x*)\f{\\x - e\\)dx, X* G K 

<mi||0||^sup||Vp(x)|| 

< midi\\e\\l sup [^•'(llxlld)! 

x£V„ 

. , |,„||,_i 11^ QJMld) \\x\\d\Q'{\\x\\d)\ ^ ,„. 

< midi||6'||^ sup f—f- sup sup x p[e) 

x(iVl' \\x\\d x(^vi' Q\W\\d) xdVl' QK\\x\\d) 



< -^ max ( (^^— j , {di(l + u)y' 1 max(|ti|, jtal) x p{9) 

for < e < 1. Therefore we have h < Cip{9), where Ci = df (1 + z^)"mi{df/(l - 
i^)}max({(l - iy)/diy^,{di{l + z/)}*^) max(|ti|, ItsD- 

Now we consider the integral outside of VJ/. We only consider \\9\\d > niax((iiz^~"'^ro,ri^ 
Then for x G V^p 

||x-6'|| > iy\\9\\ > iyd^'^\\9\\d > ro. 

Therefore we have, for < a < s 

\\xrf{\\x-9\\)dx< [ {\\x-9\\ + 11911}" f{\\x-9\\)dx 
v,9 Jv,9 



< (l + l/z^)° / \\x-9\\'^f{\\x 
Jv,9 



2X 



00 



< (1 + l/u)''cpL / r-'+'^-'dr 

Jiy\\e\\ 

= (1 + l/i/)"cpL^^^^^ 

s — a 

< C2ia)\\9\\r' , (2.11) 

where 6*2(0) = (1 + l/h')°'CpL{di/h'y^°'{s — a)~^. Hence for the second term I2, if s > 1 
and < e < 1, then I2 < C2{a)p{9) for a = 0, 1. 

We have seen that Ii and I2 are bounded from above assuming only s > 1. 

The third term I3 of 1)2. lUj) is more problematic. Write 

^3 = ii^iir / ji^ri/'(^)i/(ik-^ii)d^ 

< IMd'^if +, 

Vvcn{\\x\u<r^} Jvcn{r^<\\x\u<\\e\U} 
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+ / \\\xr\p{x)\f{\\x-9\\)dx 

= hi + h2 + h3- (say) 

We take care of 133 first. Since Q{r)r~^^ is monotone nonincreasing for r > ri, p(x)||x||^ ^ < 
/'(^)Plld ^ fo'^ W^Wd > ||^||d(> '"i)- Therefore we have, 

/33 < ii0iir""*v(ex"^*"°^ / ikr+*v(ik - o\\)dx. 

If < a + t2 < s, as in (ITTTII 

l|a;|r^*'/(lk-^||)(ix< / llxW^^^^ f{\\x-9\\)dx 
V'.c^n{||a'||^>||e|U} JVC 

<C2{a + t2)\\er/''-' 

and if a + ^2 < 0, 

Ikir^^Vdlx - e\\)dx < dv'^-^wef/'^ f f{\\x - e\\)dx 

vcn{\\x\u>\\e\U} JVC 

<d^'''-''C2{o)\\er/''-'. 

Hence I33 < C^zpiO) where C33 = d;"'^"^*^'") max(C72(a + ^2), C2(0)(ir"~*'). 
Next we consider /31. For \\6\\d > max(diz^^-'^ro, ri) and x G V^ 

f{\\x - 9\\) <L\\x- er^-' < L{u\\9\\r^-' < L{di/u\\e\Uy+^. (2.12) 

Therefore 

Jsi S WWd -t^«i V WWd / \\^\dP\x)^x- 

"'l|a:|U<»"l 

Note that by simple change of variables we have 



/ dx = Cp^dr'P"^, Cp,d = Cp n '^i ^• 

J\\x\\d<r .'I 



5^ ^ll-IU<' ,=1 

Then 



Therefore 



/ ||x||S |p(x)|dx = Cp,rf r rP+°-i|^(r)!dr 

7||a:|L<ri Jo 



Jo 



where C* = Ld^ ^v p ^Cp^d- On the other hand for 116*11^ > ri, p{9) = ^(||0|ld) is bounded 
from below as 



Q{ri)r^nen<Qm\d). 
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Therefore 

/31 < ll^lir""''"'"*^ X C*4^ r rP+--'\Q{r)\dr x p{e). 

g{ri) Jo 

Hence if s > —ti — a — p, then we can choose e > (say e= {a + p + s + ti)/4) such that 
e — a — p — s — ti <0 and hence /31 < C^ip{6) where 



Qin) Jo 



'gin) Jo 

Finahy we consider 132. Note Q{r) < £'(||6'||d)||6l||/V*i for n < r < \\e\\d and (|2?T2|) 
Then 

I32 < ||^|ir""*^L<+V-P-^||0||/-^^(||0|U) / Wxt'+'^dx 

•'ri<\\x\\a<\\9\\d 



d 

ri 



< ll^lir""''"'"*' xCj rP+'^+ti-^dr X p{e). 



Consider the integral Q = /j' "''■ rP+'^+^i ^dr. If p + a + ti < 0, then 

Q<r{^+-+P/^_t,-a-p). 

Therefore as in the case of /31, if s > — ti — a — p, then we can choose e > (say 
e = {a+p + s + ii)/4) such that e — a — p — s — ti <0 and hence 

I32 < rr'CMe) < c^2P{e) 

where C32 = {-l/{p + a + ti)}C^. Ifp + a + ti > 0, 

|P+a+ti+€3 -eg 



t\w\\d WHW 

Q = / ^P+'^+tl+^3-l-.3^^ < \U1± 

Jr, P + 



In p + a + ti + es 

for any €3 > 0. Hence 

p + a + ti + es 
If s > 1, we can choose e and £3 (say e = €3 = 1/4) such that e + £3 — s < and hence 

h2 < CS2P{0), 

where C32 = {p + a + ti + £3)^^(7^,. 

We have now confirmed that if s > max(l, —ti — a — p, + ^2); there exist e > and C = 
Ci + C'2 + C'3i + C'32 + C'33, suchthat (fT^ folds for \\e\\d > max((iiz^-Vo,ri, di(l -z^)"Vi) 
which equals to 2di max(ro,ri) for z/ = 1/2. D 
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In the next section, we need asymptotic behavior of the expectation of p{X) x KJ[X) 
where hi{e) = Hi{\\e\\d) given by ((T^ and 7 > 0. 

Corollary 2.1. Assume Fl, Bl and B2. For a = or 1, i > 1, 7 > and j = 1, . . . ,p, 

if s > max(l,7 — ti — a — p,t2 + a) and Jq rP'^'^^^\Q{r)\dr < 00, there exists e > (say 
e = min(l, s + ti + a + p — 7)/4j such that 



\r 



Ee[xfp{x)h]{x)] - 0^p{e)h]{e) < cp{e)h]{0) (2.13) 

for \\0\\(i > 2di max(ro,ri,r7o)- Moreover C does not depend on i. 
Proof. Since we have 

r^{Q{ri)H]{ri)}'/{e{r^)H]{r,)} = rje'{r,)/ Q^r^) + ^riH'M/HM, 

under Assumption B2 and by 5 of Theorem 12. 11 for ei = (s + ti +a+p — 7)/16(> 0) there 
exists TjQ such that 

ti - 7 - 7ei < r^{Q{ri)H] {!-,)}' /{e{r^)Hr<{r,)} < t^ 

for all r] > inax{rjQ,ri). Then for e = min(l, s + ti + a+p — 7)/4(> 0), (|2.13|) follows from 

Theorem O 

For Ai = max(r/o,ri), 

\p+'^-^\H^(r)Q{r)\dr < ^^./. .^ . / \p+''~'^\Q{r)\dr 



HUXi)QiXi)Jo i.w.v/i - HJ(^X,)0{X^) Jo 

which implies that C does not depend on i. D 

3. Admissibility 

In this section, we give a sufficient condition for admissibility of the generalized Bayes 
estimator with respect to a elliptically symmetric prior density g{9) = G{\\6\\d)- The 
assumptions on G are the following. 

Gl G{rj) is continuously differentiable in r/ > 0. There exist ri > 0, ti < t2 such that 

h < r]G'{r])/G{r]) < t2 for all r/ > ri. 
Gl' By redefining ri if necessary in Gl, ti and ^2 can be taken as same sign. 
Gl" G'{rj) is continuously differentiable in r] > 0. there exist ts < t4 such that t^, < 

r]G"{r])/G'{r]) < t^ for ah r? > ri. 
G2 J^ r]'P~^G{r])drj = co and there exists < 7 < 2 such that J^ r]'P~^G{r])Hl(r])drj < 00. 

G3 limr,^ovG'{v)/G{v) = to(> 1 -p)- 

FGl J^rP-^f{r)G{r)dr < 00 and J^ rP-^F{r)G{r)dr < 00. 
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We discuss some implications of these assumptions. By the assumption G3, 

/ r]P~^G{r])dr] < oo and / riP^^\G'{ri)\dr] < oo. 
Jo Jo 

From the former integrabihty and G2, the improperness of g occurs only at infinity. By 

2 of Theorem 12. II and the assumption G2, g(6)H^ {\\9\\ii) for any fixed i is integrable and 

hence becomes a proper probability density by standardization. Since Hi(-) approaches 

1 as i — > oo, g{6)Hj{\\9\\fi) is a sequence of proper densities approaching g{0), which is 

essential for using Blyth' method. 

By Gl and Lemma O G{ri) = 0{r]^^). Therefore if t2 < -p, then g{e) = G{\\e\\d) is 

a proper prior. Since we are considering an improper g{0), we assume ^2 > —p from now 

on. Moreover r]^^ /G{r]) = 0(1) by the assumption Gl. Since Hi £ RV_i, r]P~^rj^^Hf{r]) 

for ti > 2 — p is not integrable at infinity and r]'P~^G{r])Hi{r]) is not so either. Hence we 

also assume ti <2 — p. We now discuss when the integrabihty of J^ Tff^'^G{ri)H^{ri)drj in 

G2 holds and the relationship with Gl. When we take l3{rj) as in (|2.2|) . we easily see that 

there exists Li such that Hi{rj) < Li(r/ + c)~^{log(r/ + c)}~^. If t2 < 2 — p, there clearly 

exists L2 such that G{r]) < L2rf'~^ for rj > 1. Hence 

/oo /"OO 

rf^^G{r])Hl{ri)drj < L^Ls / (r/ + c)-\\o^{r] + c)}~'^dr] 

= LlL2{logil + c)}-' <oo, 

which shows that there exists 7(= 2) for the integrabihty in G2. On the other hand, if 
t2 > 2 — p then the integrabihty in G2 may not be apparent. But the condition t2 <2—p 
is not a necessary condition for the integrabihty. 

If G{r]) is regularly varying, then for any e > 0, we can choose ri, ti,t2, ^3,^4, such 
that t2 — ti < e and t4 — ts < e for all rj > ri. However in Gl we are allowing the 
case that liminf^^oo f/G'(ry)/G(?7) is strictly less than limsup^^^ r]G'{r])/G{r]). Hence we 
are dealing with a broader class of G{rj) than the class of regularly varying functions. It 
should also be noted that ^4 and ts are not always smaller than ^2 and ti, respectively. See 
Gleluk and de HaarJ (|l987l ) for the detail. 



The generalized Bayes estimator 6g with respect to the improper density g{9) is written 
as 

Jnp0f{\\x-9\Me)de 



6g{x) 



X 



In.fi\\x-e\\)gie)de 

In,{e-x)f{\\x-e\\)g{e)de 



In,f{\\x-9\\)g{e)de 
jj,,F{\\x-0\\)s/g{e)de 



In.fi\\x-e\\)gi9)d9 ' ^- ' 

which is well-defined if both J^p F{\\x — 9\\)\/g{6)d6 and Jj^p f{\\x — 0\\)g{9)d9 are integrable 
for all X. These are guaranteed by the assumption FGl and Lemma [3. II in the below. Write 

m{i;\x)= f i;{9)f{\\9-x\\)d9 
Jrp 



Y.Maruyama and A. Takemura/ Admissibility for spherical family 14 

M{iij\x) = -^ / ^{e)F{\\e - x\\)de 

where Cf = {t^p/"^ /T{p/2 + l)} j^ zP+^f{z)dz. Notice that F{-)/Cf is a probabihty density 
function because 

y-9\\-F{\\y-0\\)dy= f ||y|r{ / sfis)ds}dy 

RP Jrp J\\y\\ 

Cp / rP-i+" / sf{s)dsdr 

Jo Jr 

/•OO /"OO 

Cp / rP+i+" / tf{rt)dtdr 

/OO /"OO 

t{ / r'P+^+'^f{rt)dr}dt 

/OO /"OO 

M(V5|a;) 



Then 6g is written as 



(53(x) = x + Cf 



m{g\x) 
Note that by Gl the j'-th element of Vg is given by 

Vjg(O) - djt)j—— . 

\r\\d 

We also write 

hi{x) = Hi{\\x\\d). 

Now we state the following lemma in preparation of our main theorem. 
Lemma 3.1. 1. Assume Gl, G3 and Fl. Then 

m{g\x) ^ g{x) z/s > max(l, -ti - p, ^2) (3.2) 

mlghjlx) ^ g{x)h'J{x) i/s > max(l,7 — ti — p, ^2) (3-3) 

M{gh]\x)^g{x)h]{x) i/s > 2 + max(l,7 - ti - p, ^2) (3.4) 

M{g\\9\\^^\x) ^ gix)\\x\\f i/s>2 + max(l,l-ti-p,t2-l) (3.5) 

Migh]\\e\\2'\x) « g{x)h]{x)\\x\\f ifs > 2 + max(l,7 + l-ti-p,t2- 1). 

(3.6) 

2. Assume GV, Gl", G3 and Fl. Then 

M{Vjg\x) ^ Vjg{x) if s > 2 + max(l, -t^ - p, t4) (3.7) 

M{Vjgh'J\x)RiVjg{x)h]{x) if s > 2 + max{l,j - ts - p,U). (3.8) 
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3. M{'Vg\x)/m{g\x) is bounded in x if s > 2 + max(l, 1 — ti — p,t2 — !)■ 

Proof. When we consider the asymptotic behavior of M(-|x), that is, the expectation 
under the probabihty density F{\\6 — x\\)/Cf, we have only to substitute s for s — 2 in 
order to have corresponding results for Theorem 12.21 and Corollary 12.11 because under the 
assumption Fl there exists Li such that r'P^^~'^F{r) < Li for all r > vq. We easily see 
that (1221), ^M and (ETj) follow from Theorem O and that (ESj), (|S31), ^M and (HUni 
follow from Corollary 12.11 

Note that by the assumptions Gl and G3 there exists L2 such that \r]G'{ri)/G{rj)\ < L2 
for all rj>0. Then M{Vg\x)/m{g\x) < diM{\G'\ \x)/m{g\x) < L2diM{G/\\e\\ \x)/m{g\x). 
The value at x = is clearly bounded under the assumption FGl and the value at 
\\x\\d — > 00 is bounded by (|.''{.2|) and (|I15J) if s > 2 + max(l, 1 — ti — p,t2 — I). D 

By part 3 of Lemma l3. 11 M{'Vg\x)/m(g\x) is bounded in x and hence the risk function 
of 6g is finite because 

R{e, 6g) = E[\\X - 9 + CfM{Vg\X)/m{g\X)\\l] 

< Qm..E[\\X - e + CfM{Vg\X)/m{g\X)f] 

< 2Qmax {E[\\X - ef] + C}E[\\M{Vg\X)/m{g\X)f]}, 

where Qmax is the largest eigenvalue of Q. 

Now we state the main theorem of this paper. 

Theorem 3.1. 1. Assume Gl, G2, G3, Fl, FGl. Then the generalized Bayes esti- 
mator with respect to g is admissible if s > 2 + max(l, 1 + 7 — ti — p) and t2<2 — p. 
2. We also assume Gl ' and Gl ". Then the generalized Bayes estimator with respect 
to g is admissible if s > 2 + max(l, ?> — ti — p,t2,2 — t^ — p, ^4) and t2 > 2 — p. 

Although the moment conditions for s in the theorem above looks complicated, it is 
just from the assumptions Gl, Gl' and Gl" which make our class of G broader than 
the class of regularly varying functions. We see that the condition reduces to s > 3 for 
regularly varying functions G. Before giving a proof of the main theorem we present it as 
a corollary. 

Corollary 3.1. Suppose that G{r]) is regularly varying with index k for —p < k < 2 — p. 
Assume G3, Fl with s > 3 and FGl. 

1. Assume —p < k < 2—p andG{r]) is continuously differentiable. Then the generalized 
Bayes estimator with respect to g is admissible. 

2. Assume k = 2 — p and G{ri) satisfies 

G{v) < V^-ni^ /3(r)dr}2{,?/3(r?)}-i for r? > 1, (3.9) 

and is twice continuously differentiable. Then the generalized Bayes estimator with 
respect to g is admissible. 
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Proof. When —p < k < 2 — p, we can take ti = k — e/3 and t2 = k + e/3 for any 

< e < {3/2}(2 -p-k)\n Gl. Let 7 = /c + p + 2e/3. Clearly < 7 < 2 and it satisfies 
the integrability in G2. Since ^2 < 2 — p, we have only to apply part 1 of Theorem 13.11 

1 + 7 — ti— p = l + e and hence the moment condition for A;<2 — piss>3. 

Next we consider the case k = 2 — p. Note that there exists Li such that Hi{rj) < 
^iPiv)/!^ /3(r)dr by part 2 of Lemma 2.1. So if G(ry) satisfies (|3.9|) . we have 

r]P-^G{r])Hf{7j)dr] < Li f3{r])d7] < 00, 

which shows that the integrability in the assumption G2 is guaranteed. We can take 
ti = 2 — p — e, t2 = 2 — p + e, t3 = l— p — e and ^4 = 1 — p + e for any e > in the 
assumptions Gl, Gl". Since t2 > 2 — p, we have only to apply part 2 of Theorem 13.11 
Since 

3-ti-p = l + e, t2<l + e, 2-t3-p = l + e, U<e 

the moment condition for k = 2— p[ss>3. D 

In particular the boundary case in 1)3. 9|) by taking /5(??) in (|2.2|) was the motivating one 
for this paper. 

Corollary 3.2. Assume Fl with s > 3 and FGl. Then the generalized Bayes estimator 
with respect to H^H^ '^ 11^=0 -'^^Si (11^ IU + '^)' where n is a nonnegative integer andLog„(c) > 
0, is admissible. 

In the normal case. iBrownl ( 197ll )'s sufficient conditions for admissibility and inadmis- 



sibility are known. He showed that the generalized Bayes estimator with respect to g is 
admissible if 

/ M\l~'^^{'m[g\x)}-^dx (3.10) 

"'lklU>i 

diverges and inadmissible if (|3.10() converges. By Lemma l3. 11 we see that 

{l/2)G{\\x\U) <m{g\x) <2G{\\x\U) 

for sufficiently large \\x\\d and hence that G(r/) < rf'^ Y['^^Q'Log.i{r] + c) leads to admissi- 
bility and G{rj) > r/^~^nr=^o Logj(r/ + c)Log^(r/ + c) leads to inadmissibility. Therefore our 
sufficient condition in Theorem 13. II is very close to being necessary. 

We also notice that the prior density suggested in Corollarv l3 . 2l becomes |^| HiLo Logj(|^| + 
c) and HiLoI^oSidl^lU + '^) f or p = 12 res pectively, which are thicker than the Lebesgue 
measure. In the normal case. iBrownl ( 197lh has already pointed it out. 



Furthermore we in d icate that our moment condition s > 3 is very tight because, as 
pointed out in IPernd (1970|), admissibility requires the existence of moment one degree 



higher than what is needed for finite risk in various estimation problems. 
Now we give a proof of Theorem 13.11 
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Proof of Theorem \'J.l\ Let 5gi denote the Bayes estimator with respect to the proper prior 
density g{6)h]{6). Then the Bayes risk difference of 5g and 6gi with respect to the density 
g(9)h](9) is written as 



A 



[Rie,5g)-Ri6,6gi)]gi9)h]{0)de 



RP 



RP JRP 



\\5g,-e\\l]f{\\x-e\\)g{9)h]{e)dedx 



RP 



\X ||2 llx Il2 



\\S,^rQ] fi\\x-9\\)g{9)h]i9)d9 



RP 



^2 



^2 



2iSg - 5gi)Q' f 9fi\\x - 9\\)g{9)h]i9)d9]dx 

JRP } 

f{\\x-9\\)g{9)hU0)d9^dx 
M{Vg\x) M{V{gh2]\x) ' 



llr r ||2 

RP ^ iJRP 



RP 



RP 



m{g\x) 



m 



ighllx) 



m{ghj\x)dx 



M(yg\x) M{Vgh]\x) M{gVh]\x) 



m{g\x) m{ghj\x) m{gh'J\x 

In the same way as in lBrown and Hwand (|l982r ). we have 

M{Vg\x) M{Vgh]\x ^ 



m{gh'J\x)dx. 



Q 



A < 2C?Q 



y^;max 



RP 



+ 2C fQmax 

JRP 



m{g\x) 

M{gVh]\x 



m{gh]\x) 
= 2C}QrnUBi + Ai). (say) 

Using the Cauchy-Schwartz inequality for Aj, we have 



m{ghl\x) 

2 

m{gh]\x)dx 



m{gh]\x)dx 



7 



<Y 



RP 



M{gh]'^Vhi\x) {m{gh]\x)}^^dx 
M{gh]\x) 



RP m 



{9hl\ 



X 



-M{gh]~^\\Vhif\x)dx 



<l'l ^^^Mighr'\m\\'\x)dx 



RP rn 



{gh]\x) 



for < 7 < 2. The ratio M (ghj \x) / m{gh] \x) is bounded from above by M {g\x) /m{gh'l\x) 
and hence the value at x = is clearly bounded under the assumption FGl. By (|H.I-ij) and 
(|3.4() . we have 

||x||d-+oo ■m{gh]\x) 
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uniformly in i. This implies that there exists ci such that M{gh]\x)/m(gh'J\x) < ci for all 
X and for all i. Then 



Ai < 7^01 / M{ghJ'^\\Vhi\\'^\x)dx 
Jrp 



T^ci/ {l/Cf}F{\\x-9\\)dx g{e)Kr\9)\\Vhmfde 

JRP JRP 



.7-2, 



T^ci / g{e)Kr\e)\\vh, 

JRP 



-de. 



By 4 of Theorem O we have \\Vhi{e)\\ < 2dif3{\\e\\d)/ Jp\\^P{r)dr and together with 2 
of Theorem 12.11 \\\7hi{9)\\/hi{6) for all 9 is bounded from above by L3 independent of i. 
Therefore 

Ai < ^'ciLl f g{9)hl{9)d9, 

JRP 

which is bounded by the assumption G2. Furthermore ||V/ii(0)|| ^ as i ^ 00 by 3 
of Theorem 12.11 Therefore by the dominated convergence theorem Ai converges to as 
i ^^ 00. 

Next we consider Bi. M(V5'|a;) and MCVgHJlx) at x = are zero vectors because g 
and h] are function of \\9\\d- So the integrand of Bi is bounded around x = 0. For the 
asymptotic property of the integrand of Bi, we need to distinguish two cases: t2 < 2 — p 
and t2 > 2 — p. When t2 < 2 — p, we can bound the norm in the integrand of Bi from 
above somewhat roughly. Using ()3.2|1 . 1)3. 3|1 . (|3.5p and p.6p in Lemma IXTl and noting that 
\r]G'{r])/G{r])\ < L2 by the assumptions Gl and G3 we have 






MiV,g\x) MiV,gh]\x) 



m{g\x) 



m{ghj\x 



M{\G'\ \x) MilC'lkJlx 



<L 



m{g\x) ' m{gh]\x) 

f M{G\\9\\-^\x) M{Gh]\\9\\-'^\ 
\ 'm{g\x) m{ghj\x) 



< c\\x\ 



for all sufficiently large \\x\\d and for all i. When t2 > 2 — p, we have to bound it from 
above more strictly. By (|3.2() . 1)3. 3() . 1)3. 7|) and 1)3. 8|1 in Lemma IXTl we have 



1 



M{Vjg\x) M{Vjgh]\x) 



m{g\x) m{gh]\x) 



G'{\\x\U) xj + 0{\\x 



il-eo^ 



G{\\x\\d)\\x\U l + 0(||x||^^«) Gi\\x\U)\\x\U l + 0(||x||7") 



< c\\x 



-1-^0 



<Hi(\\x\\d) 
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for some eg, for all sufficiently large \\x\\d and for all i. Moreover m{ghj\x) < m{g\x) and 
m{g\x) < 2G{\\x\\ii) for all sufficiently large \\x\\ii by (|3.2j) . Therefore there exist Ci, C2, 
C3 and C4 such that the integrand of Bi is less than 

rmin{C7i,C2 11x11^^+*^} 1^2 < 2 - p 

{mm{C3,C^Gi\\x\\d)Hf{\\x\U)} ifts >2-p. 

Therefore Bi converges to as i —> 00 by the dominated convergence theorem. 

Finally we conffi-m that we use (j321)-(E!ni) for t2 < 2 -p and (|SI2I)-(|S3I), (EIZJ and (HUni 
for t2 > 2 — p. Note also max(l, ^2) = 1 for t2 < 2 — p in the moment condition. D 

4. The generalized Bayes estimator with respect to the harmonic prior and 
its minimaxity 

In this section, we show that the generalized Bayes estimator with respect to the harmonic 
prior has a form simple enough to check some sufficient conditions for minimaxity under 
the quadratic loss function Li{9, 6) = \\6 — 6\\'^ given in early studies. We demonstrate that 
it is minimax for some /. 

In (|3.H1 . the generalized Bayes estimator can be also written as 

mlglx) 
For p > 3 and g{9) = ||0|p~^, we have 



m{g\x)= f{\\x-e\\)\\df~'Pdd= f{h\\)\\x-iif-Ur^ 

JRP JRP 

Jo Jo 

/•oo rn 

= Cp^y / f{rt){l + 2tcos^ + t^f^PlhP'^smP-'^^dtd^ 

Jo Jo 

' / tP-^f{rt)dt + r'^ / tf(rt)dt\ 

-fP-2^(rt)l ^ + (P - 2) / tP-^F{rt)dt + F{r)] 

= Cp{p -2) [ tP''^F{rt)dt, (4.1) 

Jo 

where r = ||x||. The fifth equality in the above equation follows from the relation 

r{l + 2t cos ^ + t^f-Pl'^ sinP-2 ^d^ = B{p/2 - 1/2, 1/2) min(t2~P, 1), 
^0 
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which is proved in Lemma |4. II in the end of this section. In the same way, we have 

CfV.^M{g\x) = -xcp{p - 2) / tP-^F{rt)dt. 

Jo 

Hence the generahzed Bayes estimator is written as 6^{X) = {l — (p^:{\\X\\)/\\X\\'^)X, where 

„2Jo 



Mr) 



J^tP-^F{rt)dt 
J^tP-^F{rt)dt' 



Some properties of the behavior of <?i*(r) are easily derived as follows. 

Theorem 4.1. 1. limr-,oo (p*{r) = {p - 2)Eo{\\X\\'^) /p. 

2. (j)i:{r) is nondecreasing in r for any f. 

3. (/),f(r)/r^ is nonincreasing in r if F{t){t'^ f(t)}~^ is nonincreasing. 

Proof. 4i^{r) can be written as Jq tP~^F{t)dt/ Jq tP~''^F{t)dt and we have 



lim 0*(r) 



/o°° tP'^F{t)dt _ P - 2 /o°° tP+^fit)dt _p-2 



f^tv^^F{t)dt p f^tv~^f{t)dt p 

The derivative of </>*(r) is calculated as 



En\\\X\ 



4>'*{r) 



{r' - t^)tP--'F{t)dt, 



{J^tP-^F{t)dtrJo 
which is nonnegative for any /. The derivative of (j)*{r)/r'^ is calculated as 

-2 



d 



*(*•«/' 



f tP-^F{rt)dt] (f tP-^F{rt)dt f tP-^f{rt)dt 
tP+^firt)dt f tP-^F{rt)dt] . 



If F{t){t'^f{t)} ^ is nonincreasing, the right-hand side of the equality above is nonpositive 
by the covariance inequality. D 

Now we consider the minimaxity of 5*. We present a brief list of known sufficient con- 
ditions for minimaxity given in previous papers, for the estimator of the form (|l.ll) with 
nonnegative and nondecreasing (j){r). 



Author 


P 


Hr)/r^ 


upper bound of 


general | 


Berger (1975) 


p>3 




2(p-2)inf,6C7i^(s)//(s) 


Brandwein (1979) 


p > 4 


\ 


2(p-2)(pi?o(||X||-^))-^ 


unimodal or / is nonincreasing | 




p>4 


\ 


2p((p + 2)i5o(||X||-^))-^ 


Ralcscu ct al. (1992) 


p = 3 


\ 


0.93{Eo{\\X\r'))-' 


F{t)/f{t) is nondecreasing 


Bock (1985) 


p > 4 


\ 


2{Eo{\\X\D)-^ 


scale mixtures of multivariate normal 




p>3 


\ 


2iEo{\\X\\-^))-^ 
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In the table, U = {t > 0\f{t) > 0} and an arrow \ means nonincreasing. It is noted 
that f{t) is nonincreasing in t if F{t)/ f{t) is nondecreasing in t and that F[t)/ f{t) is 
nondecreasing in t if / is a scale mixtures of multivariate normal. 

Combining Theorem 14.11 and the table above, we can derive a sufficient condition for 
minimaxity of (^*(X) and we state it in the following theorem for p > 4. 

Theorem 4.2. 1. Assume t~'^ F (t) / f (t) is nonincreasing. 

(a) 6^ is minimax if Eq[\\X\\'^]Eo[\\X\\-'^] < 2. 

(b) Assume also f{t) is nonincreasing. Then 5^ is minimax 

if{p'-^)E,[\\X\\'']E,[\\X\\-^]<2p\ 

(c) Assume also F(t)/f{t) is nondecreasing. Then 5* is minimax 

if{p-2)Eo[\\Xr]Eo[\\X\\-^]<2p. 

2. Assume < infsg^/ F(s)//(s) < cxd. Then 5^ is minimax 
ifEo[\\Xf]<2pmiseuF{s)/f{s). 

iBergen (|l97fj ) and iBockl (|l98fj ) gave several examples of /, checked the monotonicity 



of /(t), F{t)/f{t), and t~'^ F (t) / f (t) and calculated an upper bound of (j){r). In this paper 
we give just two examples but we believe that the estimator 6^{X) is minimax for a broad 
class of spherically symmetric distributions. 

Example 4.1. We consider /(s) = s"exp(— /?s^) for a, /? > 0. We have 



oo 



^ = y^ u''+'expiPt\l-u))du, 
which is decreasing in t. By an integration by parts, we have 

f{t) ~ 2/3 ^ 2/5 t"exp(-/3i2) 

= ^ + ^/ u'^-'eMPt'{l-u'})du 

and hence mU>oF{t)/f{t) = (2/3)"^ We also have Eo{\\X\\^) = {p/2 + a/2)/p and 
E'odl^ll"^)""'^ = ip/2 + a/2 — l)/0 . Ther efore the generalized Bayes estimator is min- 
imax if a < y fo r p > 3 by IBeraen \l97,i ) 's conditions and if a > 4 — p for p > A by 
Brandwein \l979i) 's conditions regardless of /3. Hence the estimator for p > 4 is minimax 



regardless of a and (3. 

Example 4.2. We consider f{t) = exp(-tV2) - aexp(-tV{2^}) /or < a < 1, < 
6 < 1. Note that if a < b then f is unimodal and if a > b then f is not. We easily 
see that mit>oF{t)/f{t) = 1 and that Eo[\\Xf] = p(l - a6P/2+i)/(i _ a^p/S). Because 
(1 — abP/2+^ )/(l — ahP/'^) < 2 for 0<a<l,0<6<l, the generalized Bayes estimator is 
minimax bv lBeraen \l97,i) . 
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The following lemma is stated in a more general form in 3.036 of lGradshtevn and Rvzhild 
(|2000|), but it is incorrectly stated with an errata posted on the book's web page. Maruyama 
pointed out this error and he is acknowledged in the errata for 3.036. Since a derivation 
of the formula is not easily accessible, we provide our own proof. 

Lemma 4.1. For a > —1/2 and \a\ < 1, 

f (1 + 2a cos (f + a^)~" sin^" (pdcp = B{a + 1/2, 1/2). 
Jo 

Proof. Let g{ip) = (1 + 2a cos (p + a^)"^ sin if. Then we have the derivative 

,, , „ . (acosi^ + l)(cos93 + a) 
(1 + zacosip + a^)"^ 

We see that g{f) is monotone increasing from g{0) = to g{arccos{—a)) = 1 and decreasing 
from 5((arccos(— a)) = 1 to g{7T) = 0. Therefore we have 

/ {l + 2acosip + a^)~°'sm'^°'ipdip 
Jo 

((■arccos(— a) rn \ 

Jo Jarccos{— a) / 

/•arccos(— a) 

= / (l + 2acos(^ + a^)""sin^°(^d99 

Jo 

/■arccos(a) 

+ / {l-2acosp + a'^)-°'sm^°'pdp 

Jo 

= f e{d^/dt)dt+ f s''{dp/ds)ds, (4.2) 

Jo Jo 

where t = (1 + 2a cos ip + a^)^^ sin^ tp and s = (1 — 2a cos p + a^)"^ sin^ p. Here {d(p/dt) 
and {dp/ds) are calculated as 

dp/ds = ^^^ (l -{as + A {s)f) ' ' , 
where A{t) = (1 - ty/^{l - aHy/^. Let 

Then we have h'^{t) = {2tA{t)}-'^ {2 - 2a'^t^ - 2A{tf + 2B{t)], where 

Bit) = (l - {at - ^(t)}2 - {at + A{t)Y + {a^t^ - A^{t)}'^'^ ^'^ 
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= t{l-a^), 
which imphes h{t) = t~^'^(l — t)~^". Therefore we get 

the right hand side of (g^l) = / t"'h{t)dt = B{a + 1/2, 1/2). 

Jo 

n 
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